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The weak decays of B — > Kq (1430)Z + Z — (I = h,t) are investigated in Minimal Supersymmetric 
Standard Model (MSSM) and also in Supersymmetric (SUSY) SO(10) Grand Unified Models. Neu- 
tral Higgs bosons are the point of main focus in MSSM because they make quite a large contribution 
in exclusive B — > X s l + l~ decays at large tan f3 regions of parameter space of SUSY models, as part 
of SUSY contributions is proportional to tan 3 /3. The analysis of decay rate, forward-backward 
asymmetries and lepton polarization asymmetries in B — > Kq (1430)/ + /~ show that the values of 
these physical observables are greatly modified by the effects of neutral Higgs bosons. In SUSY 
SO(10) GUT model, the new physics contribution comes from the operators which are induced by 
the neutral Higgs boson penguins and also from the operators with chirality opposite to that of the 
corresponding Standard Model (SM) operators. SUSY SO(10) effects show up only in the decay 
B — > Kq^ + where the transverse lepton polarization asymmetries deviate significantly from the 
SM value while the effects in the decay rate, forward-backward asymmetries, the longitudinal and 
normal lepton polarization asymmetries are very mild. The transverse lepton polarization asymme- 
try is almost zero in SM and in MSSM model, whereas it can reach to -0.3 in SUSY SO(10) GUT 
model which could be seen at the future colliders; hence this asymmetry observable can be used to 
discriminate between different SUSY models. 



I. INTRODUCTION 

It is generally believed that the Standard Model (SM) is one of the most successful theory of the second half of the 
last century as it has passed all the experimental tests carried out for its verifications. The only missing thing that 
is yet to be verified is the Higgs boson mass, which we hope will be measured at the CERN Large Hadron Collider 
(LHC) in next couple of years. To test the SM indirectly, rare decays induced by flavor changing neutral currents 
(FCNCs) b — ► sl + l~ have become the main focus of the studies since the CLEO measurement of the radiative decay 
b — > S7 [1]. In the SM these decays are forbidden at tree level and can only be induced via loop diagrams. Hence, such 
decays will provide useful information about the parameters of Cabbibo-Kobayashi-Maskawa (CKM) matrix Q 
elements as well as the various hadronic form factors. In literature, there have been intensive studies on the exclusive 
decays B -> P(V, A) Q, fl 3, Q, fl fl Q] both in the SM and beyond, where the notions P, V and A denote the 
pseudoscalar, vector and axial vector mesons respectively. 

Despite all the success of SM no one can say that it is the ultimate theory of nature as it has many open questions, 
such as gauge hierarchy problem, origin of masses and Yukawa couplings, etc. It is known that supersymmctry (SUSY) 
is not only one of the strongest competitor of the SM but is also the most promising candidate of new physics. One 
direct way to search for SUSY is to discover SUSY particles at high energy colliders, but unfortunately, so far no SUSY 
particles have been found. Another way is to search for its effects through indirect methods. The measurement of 
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invariant mass spectrum, forward-backward asymmetry and polarization asymmetries are the suitable tools to probe 
new physics effects. For most of the SUSY models, the SUSY contributions to an observable appear at loop level due 
to the A-parity conservation. Therefore, it has been realized for a long time that rare processes can be used as a 
good probe for the searches of SUSY, since in these processes the contributions of SUSY and SM arises at the same 
order in perturbation theory 11 1. 



In other SUSY models, Neutral Higgs Bosons (NHBs) could contribute largely to the inclusive processes B — > 
XJ + l~, because part of the SUSY contributions is proportional to the tan 3 /? [l2j]. Subsequently, the physical 
observables, like branching ratio and forward-backward asymmetry, in the large tan (3 region of parameter space 
in SUSY models can be quite different from that in the SM. Motivated by the fact, similar effects in exclusive 
B — > K(K*) l + l~ decay modes are also investigated [ill, where the analysis of decay rates, forward-backward 
asymmetries and polarization asymmetries of final state lepton indicates the significant role of NHBs. It is believed 
that the problem of neutrino oscillations can not be explained in the SM. To this purpose, the SUSY SO(10) Grand 
Unified Models (GUT) 13[ has been proposed in literature. In this model, there is a complex flavor non-diagonal 
down-type squark mass matrix element of 2nd and 3rd generations which is of the order one at the GUT scale. This 
can induce large flavor off-diagonal coupling such as the coupling of gluino to the quark and squark which belong 
to different generations. In general these couplings are complex and may contribute to the process of FCNCs. The 
above analysis of physical observables in B — > K(K*) decay is extended in SUSY SO(10) GUT model where it 

has been shown that the forward-backward asymmetries as well as the longitudinal and transverse decay widths of 
the said decays are sensitive to these NHBs effect in SUSY SO(10) GUT model which can be detected in the future 



B factories 



In this paper, we will investigate the exclusive decay Bq — > Aq(1430)ZZ ( I = fi, r), where A"q(1430) is a scalar 
meson, both in the Minimal Supersymmetric Standard Model (MSSM) as well as in the SUSY SO(10) GUT model 
[13I ]. We evaluate the branching ratios, forward-backward asymmetries, lepton polarization asymmetries with special 
emphasis on the effects of NHBs in MSSM. It is known that different source of the vector current could manifest 
themselves in different regions of phase space. For low value of momentum transfer, the photonic penguin dominates, 
while the Z penguin and W box become important towards high value of momentum transfer llj . In order to search 
the region of momentum transfer with large contributions from NHBs, the above decay in certain large tan (3 region 
of parameter space has been analyzed in SuperGravity (SUGRA) and M-thcory inspired models [ijj]. We extend this 
analysis to the SUSY SO(10) GUT model ll|, where there are some primed counterparts of the usual SM operators. 
For instance, the counterparts of usual operators in B — ► X s 7 decay are suppressed by m s /mb and consequently 
negligible in the SM because they have opposite chiralitics. These operators are also suppressed in Minimal Flavor 
Violating (MFV) models QQ, however, in SUSY SO(10) GUT model their effects can be significant. The reason 



is that the flavor non-diagonal squark mass matrix elements are the free parameters, some of which have significant 
effects in rare decays of B mesons [3] . 

The main job of investigating the semi-leptonic B meson decay is to properly evaluate the hadronic matrix elements 
for B — > A'q(1430), namely the transition form factors, which are governed by the non-perturbative QCD dynamics. 
Several methods exist in the literature to deal with this problem, among which the QCD sum rules approach (QCDSR) 
{19! I20I ] is a fully rclativistic approach and well rooted in quantum field theory. However, short distance expansion 
fails in non-perturbative condensate when applying the three-point sum rules to the computations of form factors 
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in the large momentum transfer or large mass limit of heavy meson decays. As a marriage of standard QCDSR 

cure the problem 
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technique and theory of hard exclusive process, the light cone QCD sum rules (LCSR) 
of QCDSR applying to the large momentum transfer by performing the operator product expansion (OPE) in terms 
of twist of the relevant operators rather than their dimension 24j . Therefore, the principal discrepancy between 
QCDSR and LCSR consists in that non-perturbative vacuum condensates representing the long-distance quark and 
gluon interactions in the short-distance expansion arc substituted by the light cone distribution amplitudes (LCD As) 
describing the distribution of longitudinal momentum carried by the valence quarks of hadronic bound system in the 
expansion of transverse-distance between partons in the infinite momentum frame. An important advantage of LCSR 
is that it allows a systematic inclusion of both hard scattering effects and the soft contributions. Phenomenologically, 

27, Q, radiative 



LCSR has been widely applied to investigate the semi-leptonic decays of heavy hadrons 
hadronic decays 
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3l| and non-leptonic two body decays of B meson 
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In our numerical analysis for Bq — > A'q (1430) decays, we shall use the results of the form factors calculated by LCSR 
approach in Ref. 3J1, and the values of the relevant Wilson coefficient for MSSM and SUSY SO(10) GUT models 

nn 

are borrowed from Ref. [11|, uM ■ The effects of SUSY contributions to the decay rate and lepton polarization are 
also explored in this work. Our results show that the decay rates are quite sensitive to the NHBs contribution. The 
forward-backward asymmetry is zero in the SM for these decays because of the absence of the scalar type coupling, 
therefore any nonzero value of the forward-backward asymmetry will give us indication of the new physics. It is known 
that the hadronic uncertainties associated with the form factors and other input parameters have negligible effects on 
the lepton polarization asymmetries, therefore we have also studied these asymmetries in the SUSY models mentioned 
above and found that the effects of NHBs are quite significant in some regions of parameter space of SUSY. 

The paper is organized as follows. In Sec. II, we present the effective Hamiltonian for the semilcptonic decay 
B — > KqI + 1~ . Section III contains the paramcterizations and numbers of the form factors for the said decay using the 
LCSR approach. In Sec. IV we present the basic formulas of physical obscrvablcs like decay rates, forward-backward 
asymmetries and polarization asymmetries of lepton in the above mentioned decay. Section V is devoted to the 
numerical analysis of these observables and the brief summary and concluding remarks are given in Sec. VI. 



II. EFFECTIVE HAMILTONIAN 



By integrating out the heavy degrees of freedom in the full theory, the general effective Hamiltonian for b — > sZ + T 
in SUSY SO(10) GUT model, can be written as Q 

2 10 



v H=l i=3 

8 i 

+E( c q.^)q,(m) + ^ i ( / i)q; ; (/x)) , 



(i) 



where Oi(fi) (i = 1, 
the energy scale /i 



10) are the four-quark operators and Ci(/j,) are the corresponding Wilson coefficients at 



37j . Using renormalization group equations to resum the QCD corrections, Wilson coefficients 
are evaluated at the energy scale \x = m b . The theoretical uncertainties associated with the renormalization scale 
can be substantially reduced when the next-to-leading- logarithm corrections are included. The new operators Qi{p) 
(i = 1, . . . , 8) come from the NHBs exchange diagrams, whose manifest forms and corresponding Wilson coefficients 
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can be found in 38l.l39l|. The primed operators are the counterparts of the unprimcd operators, which can be obtained 
by nipping the chiralities in the corresponding unprimcd operators. It is believed that the effects of the counterparts of 
usual chromo-magnetic and electromagnetic dipole moment operators as well as semileptonic operators with opposite 
chirality are suppressed by m s /mi, in the SM, but in SUSY SO(10) GUTs their effect can be significant, since #23 
can be as large as 0.5 [3Ql- Apart from this, S^ 11 can induce new operators as the counterparts of usual scalar 
operators in SUSY models due to NHB penguins with gluino-down type squark propagator in the loop. It is worth 
mentio ning that these primed operators will appear only in SUSY SO(10) GUT model and are absent in SM and 



MSSM 



The explicit expressions of the operators responsible for B — > Kq ( 1430)Z + Z transition are given by 



07 = 


1D7T 


0' 7 


.2 


o 9 = 




o' 9 - 


l07T Z 


O10 = 


167H 


o[ 


= ^(5^P R b)(h^ 5 l) 


Qi = 




Q'i 




Q2 = 


^(SP R b)(l l5 l), 


Q'2 





with Pl 1 



derived as 



12|: 



(2) 



1 ± 75) /2. In terms of the above Hamiltonian, the free quark decay amplitude for b — > ,s l + l can be 



M(b -> sl+T) 



Gf®-, 



V2it 



VtbV t * s \cl n {s^P L b)(lYl) + C 10 {s 7ll P L b)(h^ 5 l) 
2m b C e 7 n (si<j^P R b)(iri) + C Ql {sP R b){Tl) + C Q2 (sP R b)(h 5 l) + (Ci(m&) ~ Cj(m t )) \, 



where s = q 2 and q is the momentum transfer. The operator Oiq can not be induced by the insertion of four- 
quark operators because of the absence of the Z boson in the effective theory. Therefore, the Wilson coefficient C10 
does not renormalizc under QCD corrections and hence it is independent on the energy scale. In addition to this, 
the above quark level decay amplitude can receive contributions from the matrix element of four-quark operators, 
Yln=i(l + l~ s \Oi\b) , which are usually absorbed into the effective Wilson coefficient Cg^ (fi), that one can decompose 
into the following three parts 
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C 9 e// ( M ) = C 9 ( f i) + Y SD (z, s') + Y LD (z, a% 

where the parameters z and s' are defined as z = m c /Tnb, s' = q 2 /rn 2 . Ysd{z,s') describes the short-distance 
contributions from four-quark operators far away from the cc resonance regions, which can be calculated reliably in 
the perturbative theory. The long-distance contributions Yld{z, s') from four-quark operators near the cc resonance 
cannot be calculated from first principles of QCD and arc usually parameterized in the form of a phenomenological 
Brcit-Wigncr formula making use of the vacuum saturation approximation and quark-hadron duality. Wc will neglect 
the long-distance contributions in this work because of the absence of experimental data onB^ J/ipK^ (1430). The 
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manifest expressions for Y$d{z, s') can be written as 



Y SD (z,s') = / 1 (2 jS , )(3C 1 (/x) + C 2 (M) + 3C 3 (M)+C4(/i)+3C s (/i)+C 6 (/ i )) 
s')(4C 3 (n) + 4C 4 (fi) + 3C 5 (m) + Cain)) 

-h(0, S ')(C 3 (/i) + 3C 4 (m)) + jj(3C 3 (M) + C 4 (m) + 3C 5 ( M ) + C 6 ( M )), 



(4) 



with 
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27 9 



2 arctan ■ 



in for a; = 4z 2 /s' < 1 
for x = 4z 2 /s' > 1 



/i(0,s) = rr-rln -Ins + -itt 

27 9 /i 9 9 



(5) 



Apart from this, the non-factorizable effects 0, Q, 01 from the charm loop can bring about further corrections 
to the radiative b — > transition, which can be absorbed into the effective Wilson coefficient C^* . Specifically, the 
Wilson coefficient Cff* is given by 51 1 



C 7 e// (/i)=C 7 ( M )+C 6 ^ S7 (M), 



with the absorptive part for the b — > sec — > 57 rescattering 

'2 



77 14 / 23 (Gi(a; t ) - 0.1687) - 0.03C 2 (/x) 



G 1 {x) 



x(x 2 



5a: -2) 



3x 2 ln 2 x 



(6) 
(7) 



8(x-l) 3 4(a;-l) 4 ' 

where 7/ = a s (miy)/a; s (ju), x t = m 2 /m 2 v . Here we have dropped out the tiny contributions proportional to CKM 
sector V ub V* s . In addition, C 7 eff (fi) and C 9 eff (n) can be obtained by replacing the unprimcd Wilson coefficients with 
the corresponding prime ones in the above formulae. 



III. PARAMETERIZATIONS OF MATRIX ELEMENTS AND FORM FACTORS IN LCSR 



With the free quark decay amplitude available, we can proceed to calculate the decay amplitudes for semi-leptonic 
decays of Bq — > Kq(1430)1 + 1~ at hadronic level, which can be obtained by sandwiching the free quark amplitudes 
between the initial and final meson states. Consequently, the following two hadronic matrix elements 



(K* (p)\s 7l , 75 b\Bq, [p + q)), (K* ( P )\s<r flv j 5 q v b\B q , (p + q)) 



(8) 



need to be computed as can be observed from Eq. (JT|) . The contributions from vector and tensor types of transitions 
vanish due to parity conservations which is the property of strong interactions. Generally, the above two matrix 
elements can be parameterized in terms of a series of form factors as 



(Ko(p)\ s lvJ5b\B q ,(p + q)) 
(K*(p)\sa^ 5 q v b\B q ,(p + q)) 



1 



m B + m K : 



(2p + q)n q 2 - {m% - m%) q^ f T {q 2 ) 



(9) 
(10) 
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At the large recoil region, these form factors satisfy the following relations 36l. |52| 



U(q 2 ) = 2 T Mq% /-(<? 2 )=0, (11) 
m B + m K - 



Mq 2 ) = - mb f + (q 2 ), (12) 

m B - m K , 



with the help of the equations of motion in the heavy quark limit 
momentum q^ on both side and making use of the equations of motion 



53j. Contracting Eqs. (|9lfT0]l with the four 



9^17^2) = (rn 2 - mi)^i^ 2 (13) 
9 m (V'i7m75V'2) = -(mi +5712)^175^2 (14) 



we have 



(K*(p)\s l5 b\B ql (p + q)) = [f + (q 2 )p-q + f-(q 2 )q 2 } (15) 

nil, + m s 

To calculate the non-perturbaive form factors, one has to rely on some nonperturbative approaches. Conside ring 



the distribution amplitudes up to twist-3, the form factors at small q 2 for Bq — > K^l + l" have been calculated in 36| 
using the LCSR. The dependence of form factors fi{q 2 ){i = +,—,T) on momentum transfer q 2 are parameterized in 
cither the single pole form 



1 - a % q 2 /m 2 



IM 2 ) = 1 ^'2 , (16) 



or the double-pole form 



mi 2 )=-, — 27 i i{0 i h 4/ 4 . (17) 

1 - a iq 2 /m 2 B + biq*/m% 



in the whole kinematical region < q 2 < (tob q — uik') 2 while non-perturbative parameters cij and hi can be fixed 
by the magnitudes of form factors corresponding to the small momentum transfer calculated in the LCSR approach. 
The results for the parameters aj, bi accounting for the q 2 dependence of form factors / + , /_ and fx arc grouped in 
Table H 

TABLE I: Numerical results for the parameters fi(0), a.i and bi involved in the double-pole fit of form factors (|17|l responsible 
for Bo — > Kq(14:30)11 decay up to the twist-3 distribution amplitudes of _/fo(1430) meson. 







a, 


bi 


/+ 


u - y '-0.20 


0.86t° ;« 




f- 


0.073tg;g 


2.50«;f 7 


1 09+0. 69 


h 




0.69l°j? 





IV. FORMULA FOR PHYSICAL OBSERVABLES 



In this section, we are going to perform the calculations of some interesting obscrvablcs in phenomenology like the 
decay rates, forward-backward asymmetry as well as the polarization asymmetries of final state lepton. From Eq. 
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it is straightforward to obtain the decay amplitude for Bq — > K^l + l as 

M So ^. m - = Kayo + i£aw) + !*(«)] , (is) 

where the functions and T 3 are given by 

Tl = i {<%" (%») f + (q 2 )p, + 4 "" b fe" - <%") f + (q 2 )p„ (19) 
\ / tub + nix* V / 

T 2 = i (Cio - C( ) (/ + (g 2 )p M + /-(<z 2 )g M ) - - ' , r (C Q2 C' Q2 ) (f + (q 2 )p ■ q + f-{q 2 )q 2 ) g M , (20) 

and 

T 3 = i (C Ql - C Ql ) — 1— (/ + (g 2 ) P • Q + /„(<? V) ■ (21) 

It needs to point out that the terms proportional to q^ in namely f~(q 2 ) docs not contribute to the decay 
amplitude with the help of the equation of motion for lcpton fields. Besides, one can also find that the above results 
can indeed reproduce that obtained in the SM with C[ = and T 3 = 0. 



A. The differential decay rates and forward-backward asymmetry of Bo — > K$ (1430)Z + Z 

The semi-leptonic decay Bq — > A'q (1430)Z + /~ is induced by FCNCs. The differential decay width of Bo 
(1430)Z + Z — in the rest frame of B$ meson can be written as [59f 

dT{B a -► K$(U30)l+l-) _ 1 1 



/ l^„^^(1430 );+ i-| 2 ^: (22) 



dqi (2^) 3 32m Bo J Umin 

where u = (pk*(1430) + Pi-) 2 an d (? 2 = {pi+ + Pi-) 2 ; Pir«(i430)> Pz+ andp;- are the four-momenta vectors of -R'q(1430), 
l + and respectively; |-^B _>if*(i43o)i+z- | 2 is the squared decay amplitude after integrating over the angle between 
the lepton l~ and A'q(1430) meson. The upper and lower limits of u are given by 



(^'(1430) + E l-) 2 ~ iJ E KS(1430) ~ TO K *(1430) ~ \ E l~ ~ m f-) 2 ' 



U m in - (E* K , {1430) +E*-) 2 - (yj ^(1430) ~ TO /^(1430) + \j E l~ ~ m ]-f ( 23 ) 

where the energies of ^^(1430) and l~ in the rest frame of lepton pair E x»(\43o) and are determined as 

p* . TO lo ~ m gp(i43o) -g 2 . _ g 2 (9d) 



q 2 Wq 
Collecting everything together, one can write the general expression of the differential decay rate for Bq 



K*{U30)l+r as [60 1 



dT _ G 2 F a 2 \V tb V t ** 2 



ds 3072m|7r 5 s 



l \A\ 2 (2m 2 + s) A + 12sm 2 (m| - m 2 K , - sj (CB* + C*B) + 12m 2 s 2 |C| 2 + 6t \D\ 2 (t - 4m?) 

+ |B| 2 ("(2m 2 + s) (ml - 2m B m 2 K , - 2sm^.) + (n&. - s)* + 2m 2 (m A K , + Wtm 2 K , + s 2 ) j J, (25) 
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where 

A = X(m%,m 2 K , ,s)=m% + m,. + s 2 - 2m%m 2 K , - 2m 2 K ,s - 2sm|. (26) 
The auxiliary functions arc defined as 

(CI" - C?") f + ( q 2 ) + - ^ - (C?" - C??") /t(<? 2 ) 



A = 

B = i{C w -C' 1Q )f + {q 2 ) 



c = i (c w - c[ o) f ( q 2 ) + 2roe(rot + mj) (p • */ + (s 2 ) + (c Q2 - c> Q2 ) 

D = ^^^■ qI+{q2)+q2fT{q2) ^ c ^~ c ^ (27) 

1 1 lb ~r lll s 

The forward-backward asymmetry for the decay modes Bq — ► Kq(1430)Z + Z - is exactly equal to zero in the SM 
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due to the absence of scalar-type coupling between the lepton pair, which serves as a valuable ground to test 
the SM precisely as well as bound its extensions stringently. The differential forward-backward asymmetry of final 
state leptons in different SUSY models can be written as 

dA FB (q 2 ) f 1 Q d 2 T( S , COS0) f° Q d 2 T(s, C Q S e) 



ds J Q dsd cos 9 J_ 1 dsd cos 9 

and 

)d 2 r(s,cose) rO j a d 2 r(s.cost 

> ^ '- — ' dcosP ^ 



T 1 J a dT(s,cos8) rO , 

L d cos 9 . \, — 7T- — -, «< 

* / 2\ JO dsdcos9 J — 1 asa cos t7 /'OO'V 

^fM<7 J — -jj— „d 2 r(s,cose) . f o j /,d2r(s.co S e) ' ^ ' 

L a cos , y — s-^ + , dco&9—rY — ^ 

JO dsd cos 9 J— 1 dsd cos 9 



Now putting everything together, we have 

„ ^ n/ dTa 2 G 2 F \V tb V*\ 2 X(m 2 B ,m 2 K ,, S ) 4mK, 

Am W - (V*) T^pB ^(l--J-)(AD-+A-D). (30) 

It is clear from the expressions of decay rate and forward-backward asymmetry that the contribution of the NHBs 
as well as that of the SUSY SO(10) GUT model comes in through the auxiliary functions defined in Eq. ([27|) . Hence 
these SUSY effects manifest themselves in the numerical results of these observables. 



B. Lepton Polarization asymmetries of Bo — * Kq (1430)? + ? 
In the rest frame of the lepton l~ , the unit vectors along longitudinal, normal and transversal component of the l~ 



can be defined as 



63J: 



= (0,eW) = (0, , (31) 

V Ko x p- 1 / 

,s T M = (0, e T ) = (0, e N x g L ) , 

where p- and px* are the three-momenta of the lepton l~ and ^(1430) meson respectively in the center mass (CM) 
frame of system. Lorentz transformation is used to boost the longitudinal component of the lepton polarization 
to the CM frame of the lepton pair as 
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where Ei and mi are the energy and mass of the lepton. The normal and transverse components remain unchanged 
under the Lorentz boost. 

The longitudinal (Pl), normal (P/v) and transverse (Pr) polarizations of lepton can be defined as: 

p(=F)/„\ _ ds Is c ) ds VS t __[ /oo\ 

1 () "^ = ^) + f(^ = -^ } 

where i = L, N, T and £^ is the spin direction along the leptons P. The differential decay rate for polarized lepton 
l T in Bo — > (1430)Z + P decay along any spin direction £ T is related to the unpolarized decay rate (|22|) with the 
following relation 



(34) 



We can achieve the expressions of longitudinal, normal and transverse polarizations for Bq — ► Kq(1430)1 + 1 decays 
as collected below. The longitudinal lepton polarization can be written as 

dT a 2 G 2 F \V tb V*\ 2 X 3 / 2 (m 2 , mi ,,s) 4 m? 



Similarly, the normal lepton polarization is 



and the transverse one is given by 



(m| - m 2 K , + s){A*B + AB*) - 2s{A*C + AC* 



(36) 



,2 



D , . . /d r -ia 2 G||^*| 2 AV2 (m 2 ^ S ) 4m . 
P - (S) - ^ 4096^ °—{l—f)(m%-m% + s) 



(A*D - AD*) + 2mi(B*C - BC*) 

(37) 



The appearing in the above equation is the one given in Eq. ([25)) and \(m 2 B , m 2 K ,, s) is the same as that defined 
in Eq. (261). 



V. NUMERICAL ANALYSIS 

In this section, we would like to present the numerical analysis of decay rates, forward-backward asymmetries 
and polarization asymmetries. The numerical values of Wilson coefficients and other input parameters used in our 
analysis are borrowed from Ref. Q, S| an d collected in Tables II, III and IV. In the subsequent analysis, we 

TABLE II: Values of input parameters used in our numerical analysis 



G F = 1.166 x 1(T 5 GeV" 2 |V ts | = 41.61±g;Jg x 1(T 3 
\V tb \ = 0.9991 m b = (4.68 ± 0.03) GeV 

m c (m c ) = 1.2751°;°^ GeV m s (1 GeV) = (142 ± 28) MeV 
m B o = 5.28 GeV m K , = 1.43 GeV 



will focus on the parameter space of large tan /3, where the NHBs effects are significant owing to the fact that the 
Wilson coefficients corresponding to NHBs are proportional to [pibmi/mh) tan 3 j3 (h = h°, A ). Here, one tan/3 
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TABLE III: Wilson Coefficients in SM and different SUSY models but without NHBs contributions. The primed Wilson 
coefficients corresponds to the operators which are opposite in helicities from those of the SM operators. 



Wilson Coefficients 


c e 7 ff 


C? ff 


c 9 


cS 


Cio 


Cio 


SM 


-0.313 





4.334 





-4.669 





SUSY I or II 


+0.3756 





4.7674 





-3.7354 





SUSY III 


-0.3756 





4.7674 





-3.7354 





SUSY SO(10) (A = -1000) 


-0.219 + Oi 0.039 - 0.038i 4.275 + Oi 0.011 + 0.0721? 


-4.732 - Oi 


-0.075 - 0.670i 



TABLE IV: Wilson coefficient corresponding to NHBs contributions. SUSY I corresponds to the regions where SUSY can 
destructively contribute and can change the sign of CV, but without contribution of NHBs, SUSY II refers to the region where 
tan f5 is large and the masses of the superpartners are relatively small. SUSY III corresponds to the regions where tan j3 is large 
and the masses of superpartners are relatively large. The primed Wilson coefficients are for the primed operators in eq.Q from 
NHBs contribution in SUSY SO (10) GUT model. The values in the parentheses are for the r case. 



Wilson Coefficients 


Cqi 








SM 














SUSY I 














SUSY II 


6.5(16.5) 





-6.5 (-16.5) 





SUSY III 


1.2(4.5) 





-1.2 (-4.5) 





SUSY SO(10) (A = -1000) 


0.106 + 0i 
(1.775 + 0.002i) 


-0.247 + 0.242i 
(-4.148 + 4.074i) 


-0.107 + 0i 
(-1.797 - 0.002i) 


-0.250 + 0.246i 
(-4.202 + 4. 128i) 



comes from the chargino-up-type squark loop and tan 2 /3 comes from the exchange of the NHBs. At large value of 
tan/3 the Cq 1 . compete with c\ ' and can overwhelm c\ ^ in some region as can be seen from the Tables III and IV 



12j . SUSY I corresponds to the regions where SUSY can destructively contribute and can change the sign of CV, but 



without contribution of NHBs, SUSY II refers to the region where tan/3 is large and the masses of the superpartners 
are relatively small. SUSY III corresponds to the regions where tan (3 is large and the masses of superpartners are 
relatively large. The primed Wilson coefficients are for the primed operators in Eq.([2]) from NHBs contribution in 
SUSY SO(10) GUT model. As the NHBs are proportional to the lepton mass, the values shown in the table are for 
\x case and r case (the values in parentheses of table IV). Apart from the large tan/3 limit, the other two conditions 
responsible for the large contributions from NHBs are: (i) the mass values of the lighter chargino and lighter stop 
should not be too large; (ii) the mass splitting of charginos and stops should be large, which also indicate large mixing 
between stop sector and chargino sector ll|. Once these conditions are satisfied, the process B — > X s -f will not only 



impose constraints on CV but it also puts very stringent constraint on the possible new physics. It is well known that 
the SUSY contribution is sensitive to the sign of the Higgs mass term and SUSY contributes destructively when the 
sign of this term becomes minus. It is pointed out in literature ll| that there exist considerable regions of SUSY 
parameter space in which NHBs can largely contribute to the process b — > sl + l~ due to change of the sign of CV 
from positive to negative, while the constraint on b — > is respected. Also, when the masses of SUSY particles are 
relatively large, say about 450 GeV, there exist significant regions in the parameter space of SUSY models in which 
NHBs could contribute largely. However, in these cases C7 does not change its sign, because contributions of charged 
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Higgs and charginos cancel each other. We hope that these scalar mode decays of B mesons can be used to distinguish 
between these two regions of SUSY. 

The numerical results for the decay rates, forward-backward asymmetries and polarization asymmetries of the 
lepton are presented in Figs. 1-5. Fig. 1 describes the differential decay rate of B — ► K* (1430)l + l~ , from which 
one can see that the supersymmetric effects are quite distinctive from that of the SM both in the small and large 
momentum region. The reason for the increase of differential decay width in SUSY I model is the relative change in 
the sign of , while the large change in SUSY II model is due to the contribution of the NHBs. As for the SUSY 
III and SUSY SO (10) models, the value of the Wilson coefficients corresponding to NHBs is small and hence one 
expects small deviations from SM. Similar effects can also be seen for the tauon case in Fig. lb. 

In Fig. 2, the forward-backward asymmetries for B — * K$l + l~ (I = /i,r) are presented. In SM the forward- 
backward asymmetry is zero for this decay because there is no scalar operator. However in SUSY II, SUSY III 
and SUSY SO(10) model, we have the scalar operators corresponding to the NHBs, therefore we expect the nonzero 
value of the forward-backward asymmetry. This is quite clear from the Eq. (|30p where the auxiliary function D 
corresponds to the contributions from NHBs. Fig. 2a describes the forward-backward asymmetry for B — » KqH + \i~ . 
As the forward-backward asymmetry is proportional to the lepton mass, therefore for the muons case it is expected 
to be very small compared to the tauons case. Thus the maximum value of the forward-backward asymmetry is 0.05 
in SUSY II model which is hard to be observed experimentally. However, for B — > KqT + t~ the maximum value 
of forward-backward asymmetry is around 0.35 in SUSY II model. The number of events required to observe this 
asymmetry are around 10 8 or so which are accessible at large colliders like the LHCb. When the final state leptons 
are the tauon pair, the effects of SUSY III and SUSY SO(10) are still too small to be measured experimentally. 

Fig. 3(a,b) shows the dependence of longitudinal polarization asymmetry for the B — > KqI + 1~ on the square of 
momentum transfer. The value of longitudinal lepton polarization for muon is around —1 in the SM and we have a 
slight deviation on this value for SUSY I and SUSY SO(10) model. However, in SUSY II and SUSY III model the 
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FIG. 2: Forward-backward asymmetry for the B — > K£l + l (I — ^t, r) decays as functions of q 2 . The dashed-dot, dashed-double 
dot and dashed-triple dot line represents SUSY II, SUSY III and SUSY SO (10) GUT model, respectively. 
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FIG. 3: Longitudinal lepton polarization asymmetries for the B — > KqI + 1 (I = fi,r) decays as functions of q 2 . The solid, 
dashed, dashed-dot, dashed-double dot and dashed-triple dot line represents SM, SUSY I, SUSY II, SUSY III and SUSY SO(10) 
GUT model, respectively. 



value of longitudinal lepton polarization approaches to zero in the large momentum transfer region. The reason is that 
in SUSY II model we have a large value of the differential decay rate and this suppresses the value of the polarization 
in the large q 2 region. In SUSY III though the value of the decay rate is not large, relatively small contribution comes 
from the Wilson coefficients C\q. In large q 2 region, the longitudinal lepton polarization approaches to zero in 

all the models including the SM, because the factor \(m 2 B , fn 2 K* , s) goes to zero at large value of transfer momentum. 
Similar effects can be seen for the final state tauon but the value for this case is too small to measure experimentally. 

The dependence of lepton normal polarization asymmetries for B — > KqI + 1~ on the momentum transfer squared 
are presented in Fig. 4. In terms of Eq. (|36p . one can observe that this asymmetry is sensitive to the contribution of 
NHBs in almost all the supersymmetric models. Fig. 4a shows the normal lepton polarization for B — > Kq^^ . It 
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FIG. 5: Transverse lepton polarization asymmetries for the B — > KqI + 1 (I — fi, r) decays as functions of q 2 



can be seen that P/v changes its sign in the case of SUSY III model and this is due to the contribution from NHBs. 
Now for SUSY II model, though large contributions from NHBs but it is overshadowed by the opposite sign of C^ e ^ 
and Cg . As the normal lepton polarization is proportional to the lepton mass, for t + t~ channel, it is expected that 
one can distinguish between different SUSY models, which can be seen from the Fig. 4b. Again due to same reasons 
as for the muons case, the normal lepton polarization changes its sign in SUSY III model. 

Fig. 5 shows the dependence of transverse polarization asymmetries for B — ► K$l + l~ on the square of momentum 
transfer. From Eq. (|3"T]) we can see that it is proportional to the imaginary part of the Wilson coefficients which 
are negligibly small in SM as well as in SUSY I, SUSY II and SUSY III models. However, complex flavor non- 
diagonal down-type squark mass matrix elements of 2nd and 3rd generations arc of order one at GUT scale in SUSY 
SO (10) model, which induce complex couplings and Wilson coefficients. As a result, non zero transverse polarization 
asymmetries for B —> KqI + 1~ exist in this model. Now for channel, the value of transverse polarization 

asymmetry is around —0.3 in almost all value of q 2 except at the end points. Experimentally, to measure (Pt) of 
a particular decay branching ratio B at the no level, the required number of events are N = n 2 /(B(P T ) ) and if 
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(Pt) ~ 0.3, then the required number of events are almost 10 s for B decays. Since at LHC and BTeV machines, 
the expected number of bb production events is around 10 12 per year, so the measurement of transverse polarization 
asymmetries in the B — > K^l + l~ decays could discriminate the SUSY SO(10) model from the SM and other SUSY 
models. 

VI. CONCLUSION 

We have carried out the study of invariant mass spectrum, forward-backward asymmetry, polarization asymmetries 
of semileptonic decays B -> K^(U30)l + l~ (Z = fx,r) in SUSY theories including SUSY SO(10) GUT model. Partic- 
ularly, we analyzed the effects of NHBs to this process and the main outcomes of this study can be summarized as 
follows: 

• The differential decay rates deviate sizably from that of the SM especially in the large momentum transfer region. 
These effects are significant in SUSY II model where the value of the Wilson coefficients corresponding to the 
NHBs is large. However, the SUSY SO(10) effects in differential decay rate of B -> iCj(1430)/+r (I = n, t) are 
negligibly small . 

• The forward-backward asymmetry for the decay B — > K^l + l" is zero in the SM because of the missing of scalar 
operators in SM. Hence, the SUSY effects show up and the maximum value of the forward-backward asymmetry 
is around 0.35 for B — > K$t + t~ in SUSY II model. When the final state leptons are the tauon pair, the effects 
of SUSY III and SUSY SO (10) are still too small to be measured experimentally 

• The longitudinal, normal and transverse polarizations of leptons are calculated in different SUSY models. It 
is found that the SUSY effects are very promising which could be measured at future experiments and shed 
light on the new physics signal beyond the SM. The transverse polarization asymmetry is the most interesting 
observable to look for the SUSY SO(10) effects where its value is around 0.3 in almost all the q 2 region. It is 
measurable at future experiments like LHC and BTeV machines where a large number of bb pairs are expected 
to be produced. 

In short, the experimental investigation of observables, like decay rates, forward-backward asymmetry and lepton 
polarization asymmetries in B — > A'q (1430)Z + Z~ (7 = /i,r) decay will be used to search for the SUSY effects, in 
particular the NHBs effect, encoded in the MSSM as well as SUSY SO(10) models. 
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